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Abstract 

We provide an equation of state for high density supernova matter by applying a momentum- 
dependent effective interaction. We focus on the study of the equation of state of high-density 
and high-temperature nuclear matter containing leptons (electrons and neutrinos) under the 
chemical equilibrium condition. The conditions of charge neutrality and equilibrium under 
/3-decay process lead first to the evaluation of the lepton fractions and afterwards the evalu- 
ation of internal energy, pressure, entropy and in total to the equation of state of hot nuclear 
matter for various isothermal cases. Thermal effects on the properties and equation of state 
of nuclear matter are evaluated and analyzed in the framework of the proposed effective 
interaction model. Since supernova matter is characterized by a constant entropy we also 
present the thermodynamic properties for isentropic case. Special attention is dedicated to 
the study of the contribution of the components of /3-stable nuclear matter to the entropy 
per particle, a quantity of great interest for the study of structure and collapse of supernova. 

PACS number(s): 21.65.+f, 21.30.Fe, 24.10.Pa, 26.60.+C, 26.50.+X, 26.60.Kp 

Keywards: Hot Nuclear Matter, Effective Interaction, Equation of State, Nuclear Sym- 
metry Energy, Proton Fraction, Supernova, Neutron Star. 



1 Introduction 



Knowledge of the properties of the equation of state (EOS) of hot asymmetric nuclear matter is 
of fundamental importance to understand the physical mechanism of the iron core collapse of a 
massive star which produces a type-II supernova, and the rapid cooling of a new born hot neutron 
star. Additionally, the EOS defines the chemical composition, both qualitative and quantitative, 
of the hot nuclear matter. [H [21 [Sj H] . Supernova explosions and neutron stars provide a unique 
laboratory where the EOS of nuclear matter can be investigated. A great opportunity to explore 
the EOS and properties of dense neutron-rich matter is available at the accelerator facility at GSI, 
and in future, will be also available through the high-energy radioactive beams at the planned 
Facility for Rare Isotope Accelerator (FRIA) |5l |6] . 

There is a wealth of existing literature regarding the EOS of supernova matter [11 HJ [TJ El El 
[10l[ni[T2l[T3l[Tl[l5l[T6l[I71[l8l[T9l[^ Supernova 
matter which exists in a collapsing supernova core and eventually forms a hot neutron star at 
birth is another form of nuclear matter distinguished in the participation of degenerate neutrinos 
and electrons [27]. It is characterized by almost constant entropy per baryon S = 1 — 2 (in units 
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of the Boltzmann constant ks) throughout the density n and also by a high and almost constant 
lepton fraction Yi = 0.3 — 0.4 in contrast with ordinary neutron star matter where S = and 
Yi < 0.05. These characteristics are caused by the effects of neutrino-trapping which occurs in the 
dense supernova core where a neutron star is formed. 

This paper is a continuation of our previous work concerning the EOS of hot /3-stable nuclear 
matter in cases where neutrinos have left the system [33]. More specifically, in order to study 
the properties and the EOS of hot nuclear matter, a momentum-dependent effective interaction 
model (MDIM) has been applied, one which is able to reproduce the results of more microscopic 
calculations of dense matter at zero temperature and which can be extended to finite temperature 
[21 ESI EH ES] • The main incentive for the present study is the fact that only few calculations of 
the equation of state of the supernova matter at high densities are available, although at lower 
densities {n < no) (where no = 0.16 fm~^ is the saturation density) reliable results are already 
available. For our purposes, we have applied a model which, in comparison to previous models, 
has the specific property that the temperature affects not only the kinetic, but also the interaction 
part of the energy density . In this way, we are able to simultaneously study thermal effects on the 
kinetic part of the symmetry energy and symmetry free energy, in addition to the interaction part 
of the above quantities [31]. This is significant in the sense that the density dependent behavior 
of the symmetry energy and symmetry free energy strongly influence the values of the proton 
fraction and as a consequence the composition of hot /3-stable nuclear matter. 

Our focus of interest is on the study of dense supernova matter. It has been speculated that 
matter at densities up to about n = Atiq may be present in the core collapse of type-II supernova 
[1]. The present work can also be applied to the study of a neutron star at its birth, which is 
of particular interest as such a star creates a new form of matter under extreme conditions. In 
particular, proto-neutron stars are identified as a final stage of a supernova collapse. At this stage, 
a proto-neutron star is hot and composed of the so-called supernova matter. 

In addition, we examine the two findings of the previous work of Takatsuka et al. [271 ES] 
concerning supernova matter. The first one is concerned with the finding that the population 
of the components is remarkably constant both in baryon density n and temperature T and the 
proton fraction Yp is very large (e.g., Yp ~ 0.3 for Yi = 0.4) in contrast with that of neutron star 
matter. The second one concerns the finding that the EOS of dense supernova matter is by far 
stiffer than that of neutron star matter and correspondingly, hot neutron stars at birth are not 
only "fat" but hot as well compared to usual cold neutron stars. We broaden our study further 
by examining the influences of the temperature on the stiffness of EOS compared to the cold case. 

The article is organized as follows. In section II the model and relative formulas are discussed 
and analyzed. Results are reported and discussed in section III, whereas the summary of the work 
is given in section IV. 

2 The model 

The model we use here, which has already been presented and analyzed in our previous papers 
[32l [33] EH ES] , is designed to reproduce the results of the microscopic calculations of both nuclear 
and neutron-rich matter at zero temperature and can be extended to finite temperature [2]. We 
provide the main characteristics of the model as follows: 

The energy density of the asymmetric nuclear matter (ANM) is given by the relation 

e(n„. Hp, T) = efci„(n„, T) + efi„(np, T) + Vintiun, Up, T), (1) 

where (n^) is the neutron (proton) density and the total baryon density is n = n„ + Hp. The 
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contribution of the kinetic parts are 



CfcinK, T) + eli^{np, T) = 2J ^^^^3 (/„(n„, k, T) + fp{np, k, T)) , 



(2) 



where fr, (for t = n,p) is the Fermi-Dirac distribution function. 

Including the effect of finite-range forces between nucleons, the potential contribution is pa- 
rameterized as follows [2] 
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In Eq. ([3]), / is the asymmetry parameter (/ = (n„ — np)/n) and u = n/uo, with Uq denoting 
the equilibrium symmetric nuclear matter density, Hq = 0.16 fm~^. The asymmetry parameter I 
is related to the proton fraction Yp by the equation / = (1 — 2Yp). The parameters A, B, a, Ci, C2 
and B' which appear in the description of symmetric nuclear matter are determined in order that 
E[n = uq) — mc^ = —16 MeV, no = 0.16 fm~^, and the incompressibility are K = 240 MeV. The 
additional parameters xq, X3, Zi, and Z2 used to determine the properties of asymmetric nuclear 
matter are treated as parameters constrained by empirical knowledge [2]. 

The first two terms of the right-hand side of Eq. ([3]) arise from local contact nuclear interaction 
which lead to power density contributions as in the standard Skyrme equation of state. The first 
one corresponds to attractive interaction and the second one to repulsive. They are assumed 
to be temperature independent. The third term describes the effects of finite range interactions 
according to the choice of the function g{k,Ai), and is the temperature dependent part of the 
interaction [5]. This interaction is attractive and important at low momentum, but it weakens 
and disappears at very high momentum. The function, g{k, Aj), suitably chosen to simulate finite 
range effects, has the following form 




-1 



(5) 



where the finite range parameters are Ai = 1.5A;^ and A2 = 3A;^ and k^p is the Fermi momentum 
at the saturation point hq. 

The energy density of asymmetric nuclear matter at density n and temperature T, in a good 
approximation, is expressed as 



e(n, T, /) = e(n, T, / = 0) + esym{n, T, I), 



(6) 



where 



e,yUn, T, I) = nl'E%l{n, T) = {E^^'^^in, T) + ^^^^^(n, T)) . 



(7) 



In Eq. (I7j) the nuclear symmetry energy El°y^{n^T) is separated into two parts corresponding to 
the kinetic contribution E^y^{n,T) and the interaction contribution Elyl^{n,T). 
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From Eqs. and ([7]) and setting / = 1, we find that the nuclear symmetry energy T) 
is given by 

E%l{n, T) = ^ (e(n, T, J = 1) - e(n, T, 7 = 0)). (8) 

Thus, from Eq. ([8]) and by a suitable choice of the parameters Xq, x^, Zi and Z2, we can obtain 
different forms for the density dependence of the symmetry energy Elyl^{n,T). 

It is well known that the need to explore different forms for El°l^{n,T) stems from the un- 
certain behavior at high density [2]. The high-density behavior of symmetry energy is the worst 
known property of dense matter [371 EH EH], with different nuclear models giving contradictory 
predictions. Recently, the density dependence of the symmetry energy in the equation of state of 
isospin asymmetric nuclear matter has been studied using isoscaling of the fragment yields and the 
antisymmetrized molecular dynamic calculation |40]. It was observed that the experimental data 
at low densities are consistent with the form of symmetry energy, Esym{u) ~ 31.6m°'^^, in close 
agreement with those predicted by the results of variational many-body calculations. In Ref. [50] 
it was also suggested that the heavy ion studies favor a dependence of the form Esymiu) ~ 31. Gm^, 
where 7 = 0.6 — 1.05. This constrains the form of the density dependence of the symmetry energy 
at higher densities, ruling out an extremely "stiff' and "soft" dependence [40 ]. 

In a previous study conducted with isospin dependent Boltzmann-Uehling-Uhlenbeck transport 
calculations, Chen et al. |41j have shown that a stiff density dependence of the symmetry energy 
parameterized as Esym{u) ~ 31.6m^ °^ clearly explains the isospin diffusion data [42] from NSCL- 
MSU (National Superconducting Cyclotron Laboratory at Michigan State University). 

In the present work, since we are primarily interested in the study of thermal effects on the 
nuclear symmetry energy and free energy, we choose a specific form for this, enabling us to 
accurately reproduce the results of many other theoretical studies [121 El . In Ref. [43] the authors 
carried out a systematic analysis of the nuclear symmetry energy in the formalism of the relativistic 
Dirac-Brueckner-Hartree-Fock approach using the Bonn one-boson-exchange potential. In a very 
recent work, [44] the authors applied a similar method as in Ref. [43] for the microscopic predictions 
of the equation of state in asymmetric nuclear matter. In this case Esym{u) is obtained with the 
simple parametrization Egymiu) = Cu^ with 7 = 0.7—1.0 and C ~ 32 MeV. The authors concluded 
that a value of 7 close to 0.8 gives a reasonable description of their predictions although the use 
of different functions in different density regions may be best for an optimal fit [44]. The results 
of Refs. [43], m] are well reproduced by parameterizing the nuclear symmetry energy according to 
the following formula 

EX^{n,T = Q) = l?>u'/^ + nF{u), (9) 

where the first term of the right part of Eq. ^ corresponds to the contribution of the kinetic 
energy and the second term to the contribution of the interaction energy. 

For the function F{u), which parameterizes the interaction part of the symmetry energy, we 
apply the following form 

F{u)=u. (10) 

The parameters Xq, X3, Zi and Z2 are chosen so that Eq. ([8]), for T = reproduces the results of 
Eq. for the function F{u) = u. 

In general, in order to obtain different forms for the density dependence of Esym{n), the function 
F{u) can be parameterized as follows ^ 

F{u) = ^/Il, F{u)=u, F{u) = 2u^/{l + u). (11) 

It is worthwhile to point out that the above parametrization of the interacting part of the 
nuclear symmetry energy is extensively used for the study of neutron star properties [2], US] as well 
as the study of the collisions of neutron- rich heavy ions at intermediate energies [IS] HT] . For a 
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very recent review of the applications of the proposed momentum dependent effective interaction 
model and the specific parametrization of it see Ref. [18] (and references therein). 



2.1 Thermodynamic description of hot nuclear matter 

In order to study the properties of nuclear matter at finite temperature, we need to introduce the 
Helmholtz free energy F which is written as [l9l [50] 



F(n, T, /) = E{n, T, I) - TS{n, T, I). (12) 

In Eq. (1121) . E is the internal energy per particle, E = e/n, and S is the entropy per particle, 
S = s/n. From Eq. (1121) it is also concluded that for T = 0, the free energy F and the internal 
energy E coincide. 

The entropy density s has the same functional form as that of a non-interacting gas system, 
given by the equation 

5.(n,T,/) = -2 |^[/,ln/. + (l-/,)ln(l-/,)], (13) 

while the pressure and the chemical potentials defined as follows [19], |50] 



^ ^dn )^^^; ^9n) gy^^,^- ^^^^ 

At this point we shall examine the properties and the EOS of nuclear matter by considering 
an isothermal process. In this case, the pressure and the chemical potentials are related to the 
derivative of the total free energy density f = F/n. More specifically, they are defined as follows 

The pressure P can also be calculated from the equation [49l [50] 

P = Ts-e + Y.fiini. (16) 

i 

It is also possible to calculate the entropy per particle S{n, T) by differentiating the free energy 
density / with respect to the temperature, that is 

S(«.T) = -W . (17) 



^ dT 

The comparison of the two entropies, that is from Eqs. f[T3|) and f[T7|) . provides a test of the 
approximation used in the present work (see Ref. [33]). It is easy to demonstrate by applying 
Eq. f[T^ that (see for a proof [IS] as well as [ST] ) 



'dF\ _ ( dF' 
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We can define the symmetry free energy per particle Fsym{n,T) by tlie following parabolic 
approximation (see also [STl [52] ) 

F{n, T, I) = F{n, T, / = 0) + I^Fsym{n, T) = F{n, T, / = 0) + (1 - 2YpfFsy^{n, T), (19) 

where 

Fsymin, T) = F{n, T,I = 1)- F{n, T, / = 0). (20) 

It is worth noting that the above approximation is not valid from the beginning, but one needs to 
check the validity of the parabolic law in the present model before using it. In Ref. [33] we have 
proved the validity of the approximation f|T9|) . 

Now, by applying Eq. f|T9l) in Eq. f|T8|) . we obtain the key relation 

A = /i„ - = 4(1 - 2Yp)Fsym{n, T). (21) 

The above equation is similar to that obtained for cold nuclear matter by replacing Esym{n) with 
Fsym{,n,T). The finding that, for both quantities {E{n,T,I) and F{n,T, I)), the dependence of 
the asymmetry parameter / can be approximated very well by a quadratic dependence leads to 
the conclusion that the entropy S{n,T,I) must also exhibit quadratic dependence of / that is 

Sin, T, /) = Sin, T,I = 0) + I'^Ssy^n, T) (22) 

where ^ 

Ssymin,T) = Sin,T,I = 1) - Sin,T,I = 0) = —iEsymin,T) - Fsymin,T)). (23) 

In order to check the validity of the parabolic approximation (|23l) . we plot in Fig. 1 the difference 
Sin,T,I = 1) — Sin,T,I = 0) as a function of P at temperature T = 10 and T = 30 MeV for 
two baryon densities, i.e., n = 0.2 fm~^ and n = 0.4 fm~^. It is thus evident that in a good 
approximation, an almost linear relation holds between S'(n, T, / = 1) — Sin,T,I = 0) and I^, 
especially for low values of P. 

Also noteworthy in the present model is that due to temperature dependence of the interaction 
part of the energy density, the temperature affects both the kinetic part contribution on the 
entropy 5" and the interaction part. Hence, the total entropy per baryon must be written as follow 

Stot = Skin + Sint- 



2.2 /^-equilibrium in hot proto- neutron star and supernova 

Stable high density nuclear matter must be in chemical equilibrium with all type of reactions, 
including the weak interactions in which (3 decay and electron capture take place simultaneously 

n — > p + e~ + Ue, p + e~ — > n + u^- (24) 

Both types of reactions change the electron per nucleon fraction, and thus affect the equation 
of state. In a previous study, we assumed that neutrinos generated in these reactions left the 
system [33|. The absence of neutrino-trapping has a dramatic effect on the equation of state and 
is the main cause of a significant reduction in the values of the proton fraction Yp |27l 136] . 

The chemical equilibrium of reactions (^M can be expressed in terms of the chemical potentials 
for the four species 

/in + fJ-u, = yUp + l^e- (25) 

The charge neutrality condition provides us with the equation 

Yp = Fe, (26) 
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while the total fraction of leptons is given by 

Yi = Y, + y;, 

Now, according to Eqs. fl2T]) and (125|) we have 



(27) 



yWe - /i,., = /in - /ip = 4 (1 - 2Fp) Fsym{,n, T). (2. 
Moreover, the leptons (electrons, muons and neutrinos) density is given by the expression 

dk 
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(27r)= 



1 + exp 



^Jh?k'^c'^+mfc*-^^l 
T 



(29) 



where g stands for the spin degeneracy {g = 2 for electrons and muons and = 1 for neutrinos). 
One can self-consistently solve Eqs. (1271) . (125]) and fl2^ in order to calculate the proton 

fraction Yp{= Kg), the neutrino fractions Y,^^, as well as the electron chemical potential /Xg as a 
function of the baryon density n for various values of temperature T. 

Depending on the form of the symmetry energy, muons can appear at high density. Prakash 
IIS] has shown that the more rapidly F{u) increases with density, the lower the density at which 
muons appear. For example, with F{u) = u, muons appear at m ~ 4, while with the choice 
F{u) = a/u, muons cannot appear till m ~ 8. However, the presence of muons has very little effect 
on the electron lepton fractions (compared to the case without the inclusion of muons) since 
remains extremely small (~ 10~"^) over a wide range of densities [15]. Thus, we do not include the 
contribution of muons in our study. 

The next step is to calculate the energy density and pressure of leptons given by the following 
formulae 



ei{ni,T) 



^ ' ^ 3 (27r)3 
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(2vr) 



h^k'^c'^ + mfc^ dk 



1 + exp 
1 
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(30) 



(31) 



The entropy density s has the same functional form as that of a non-interacting gas system, 
given by the equation 



/f} k" 
^ [/, In /, + (l-/Oln(l -/,)]. 



(32) 



The equation of state of hot nuclear matter in /5-equilibrium (considering that it consists of 
neutrons, protons, electrons and neutrinos) can be obtained by calculating the total energy density 
etot as well as the total pressure Ptot- The total energy density is given by 



etot{n,T,I) = eb{n,T,I) + J2 e/(n,T,/), 



(33) 



where eh{n,T,I) and ei{n,T,I) are the contributions of baryons and leptons respectively. The 
total pressure is 

Ptotin, T, I) = Ptin, T,I)+ Pi{n, T, I), (34) 
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where Ph{n,T,I) is the contribution of the baryons (see Eq. fll6l) ) i.e. 

Pb{n,T,I) =T Sr{n,T,I)+ ^ nrfir{n,T, I) - eb{n,T, I), (35) 

T=p,n T=n,p 

while Pi{n, T, I) is the contribution of the leptons (see Eq. (13T!) ). From Eqs. (l33l) and (IMI) we can 
construct the isothermal curves for energy and pressure and finally derive the isothermal behavior 
of the equation of state of hot nuclear matter under /3-equilibrium. 



3 Results and Discussions 



We calculate the equation of state of hot asymmetric nuclear matter by applying a momentum 
dependent effective interaction model describing the baryons interaction. We consider that nuclear 
matter contains neutrons, protons, electrons and neutrinos under /3-equilibrium and charge neu- 
trality. The key quantities in our calculations are the proton fraction Yp and also the asymmetry 
free energy defined in Eq. (l20l) . It is worth pointing out that since the supernova explosion itself 
is a dynamic phenomenon, the chemical composition of matter changes according to the evolution 
of the star all the time |2Sj. During supernova explosion, the chemical composition of matter 
reaches equilibrium not in the whole star but locally. In our present work we assume matter in 
the chemical equilibrium for simplicity in order to analyze the properties of hot neutron star and 
supernova matter. 

The validity of the parabolic approximation f|T9l) has already been tested and presented in 
our previous work [33]. Fsym{u,T), for various values of the temperature T, was derived with a 
least-squares fit to the numerical values according to Eq. fl20l) and has the form [33] 



sym' 
Psymi^^i T 

F (u- T 



-0) 

u;T = 10) 
20) 
30) 



13^2/3 + 17m 

3.653 + 28.018m- 1.5126^2 
5.995 



0.185^3 -O.OIOm^ 



26.157n - 0.827^2 + 0.068u^ - 0.002M^ 

0.193u=^ + 0.014M^ 
0.289m=^ + 0.018m^ 



13.200 + 21. 267u + 0.800m^ 
21.087+ 17.626m + 1.645m^ 



(36) 



0, is included as well. In that case Fgym coincides with Egym- 



where the case with T 

Firstly, in order to clarify the contribution of the three terms of the potential energy density, 
we plot the terms as a function of the baryon density, in Fig. 2(a). In that figure we have that 
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The first term corresponds to attractive interaction where the second corresponds to 
repulsive interaction and dominates for high values of n (n > 0.6 fm~^). Both of the above terms 
are temperature independent. The third term V'^ contains the momentum dependent part of 
the interaction, corresponds to attractive interaction and its main contribution is to compete the 
repulsive interaction of for high values of n and as a consequence avoid acausal behavior of 
the EOS at high densities. The term V'^ consists of two finite range terms, one corresponding to 
a long-range attraction and the other to a short-range repulsion. 
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Thermal effects on the momentum dependent term V*-^ are displayed in Fig. 2(b). The con- 
tribution of V'" is plotted for various values of T. It is concluded that thermal effects are more 
pronounced for high values of T (T > 10 MeV) leading to a less attractive contribution. More 
precisely, we find that for small values of n (i.e. n = 0.15 fm~^ ) V^' increases (compared to the 
cold case T = 0) 3% — 20% for T = 10 — 30. For higher values of n the increase is even less. 

The outline of our calculations procedure approach is the following: For a fixed value of baryon 
density density n, lepton fraction Yi, temperature T and initial trial value of proton fraction Yp 
(= Ye), Eq. fl29|) is solved in order to calculate the chemical potential fie- The knowledge of fie 
allows the calculation of fii,^ from Eq. (!28|) which may then be used to infer the neutrino fraction 
Y,^^ from Eq. (129]) . Finally, a new value of proton fraction Yp (= Ye) is taken from equation 
Ye = Yi — Yy^ and the procedure is repeated from the beginning until a convergence is achieved. 

In Fig. 3 we plot the fraction of electrons Ye and neutrinos Y^,^ versus the baryon density n 
for lepton fraction Yi = 0.3 and Yi = 0.4 and for various values of T. It is concluded that thermal 
effects are important, both for electron and neutrinos fractions for low values of the baryon density 
n i.e. n < 0.4 fm~^. Ye is an increasing function of T and consequently Y^^ is a decreasing function 
of T. For higher values of n, the thermal effects on lepton's fraction are unimportant. 

At this point, following the discussion of Takatsuka et al. [27], we attempt to extend the 
discussion concerning the dependence of equilibrium fraction Ye{= Yp) on the baryon density as 
well as on the nuclear symmetry energy. We ignore the temperature effect to clarify the situation. 
Actually, the situation does not change by including finite temperature effects. The energy per 
baryon of supernova matter Esm and cold neutron star matter Ej^s are expressed as function of 
n and Yp (see also ref. [27j) as 

EsM{n,Yp) = Eb{n,Yp) + Ee{n,Yp) + E,^{n,Yp) (38) 

= Ekin, Yp = 0.5) + EsyMil - 2Fp)^ + 253.6^^/^1;^/^ + 319.516u^/\Yi - Yp)^/^ 

ENs{n,Yp) = Eh{n,Yp) + Ee{n,Yp) (39) 

= Ebin, Yp = 0.5) + Esymin)il - 21^)' + 253.6u^/%^/^ 

where the symmetry energy Esym{n) is parameterized according to Eq. f]TT]) . Esym{n) is plotted 
in Fig. 4(a) for the three different parametrizations. In the same figure we have included recent 
results provided in reference ^1] achieved by performing microscopic calculations in asymmetric 
nuclear matter. In this case Esym[n) is obtained with the simple parametrization 

Esymiu) = CU^ 

with 7 = 0.8 and C = 32 MeV. It is obvious that the results of the above parametrization, 
correspond very well with the parametrization F{u) = u which is proposed here. 

The equilibrium proton fraction Yp is calculated by solving the equation dEsM,Ns/dYp = 
for various values of the density n, Esym{n) and Yi = 0.4 for supernova matter. The results are 
presented in Fig. 4(b). In the case of cold neutron star matter, Yp depends strongly on both the 
baryon density and the values of the Egymin). This is not the case for supernova matter where 
the effect of nuclear symmetry energy in determining Yp is less important than in cold neutron 
star matter. In addition, Yp, for a fixed parametrization of F{u) is almost constant with respect 
to n. 

Fig. 5 displays thermal effects on the chemical potential of leptons for Yi = 0.3 and Yi = 0.4. 
In all cases, fit is a slightly decreasing function of T. In fact, the important quantity for our 
calculations is the difference fi = fie — fiy^ which is strictly constrained from Eq. (128]) . So, it is 
appropriate to check the validity of Eq. (]T^ at least for proton fraction Yp ^ 0.3 (or P ^ 0.16). 
We found in our previous work [53] (Fig. 1), that Eq. (1191) is accurate for the values of the proton 
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fraction which are under consideration in the present work. The use of the formula fl28l) is very 
useful since it can greatly simplify the coupled equations used for the construction of the EOS. 
We mention here that to our knowledge, the above treatment has never been applied for the study 
of the EOS of supernova matter and has been applied for the first time in the present work. 

In Fig. 6 we plot the contribution of the baryons Sb, leptons Si and the total Stot to the 
entropy per baryon. In all cases, is a decreasing function of the baryon density n. Temperature 
affects appreciably both baryon and lepton contribution. It should be noted that the contribution 
of baryons 5*;, may be written as St, = Skin + Sint, where the term Skin originates from the 
temperature effect on the kinetic part of the energy density and Sint reflects thermal effects on 
the potential energy density. More precisely, entropy density s, according to equation (fT3ll . is 
an increasing function of the diffuseness of the Fermi-Dirac distribution fr{n, k,T). As indicated 
in our previous work [33] (Fig. 11), the effect of the diffuseness of the distribution /^(n, /c,T) is 
pronounced for low values of the baryon density n and for high values of T. But as we have shown 
above, thermal effect on the momentum dependent term V'" is important for low values of n and 
also high values of T. Therefore, we conclude that the therm V'-^ has a more pronounced effect 
on the entropy density s mainly for low values of n. For higher values of n the contribution of V'" 
on s is less important. 

Furthermore, the lepton contribution Si is an increasing function of the lepton fraction Yi, 
while the baryon contribution is almost independent by YJ. For the electron and neutrino entropy 
density, our present results can be reproduced with good accuracy, at least for low values of T, by 
applying the analytical formula used by Onsi et al. [2S1 130] 

Se = Ij^T, /xe = hc{3^Y,ny/\ (40) 

s.. = l^sT, /x.^ = hc{6n'Y^ny/\ (41) 

According to the above formula, the specific contribution of the leptons (electrons and neutrinos) 
to the entropy per baryon has the form 

Se,u. = Se,Jn ~ T. (42) 

Eq. (142|) gives an nice explanation for the density and temperature dependence of Si presented in 
Fig. 6(b). 

In Fig. 7, we display the contribution to internal energy E from baryons Eb and leptons Ei 
for Yi = 0.3 and Yi = 0.4 and for various values of T. The most striking aspect is that the lepton 
energy, Ei = E^. + E,^^, dominates in the internal energy of the matter up to n ~ 0.6 fm~^ (for 
Yi = 0.3) and n ~ 0.8 fm~'^ (for Yi = 0.4). This is a characteristic of the supernova matter and is 
in remarkable contrast with the situation of cold neutron star matter [27] . The contribution from 
baryon Eb gets larger with the increase of n and is comparable with Ei for high values of n. 

It is worth pointing out that the above situation depends on the combination of the stiffness 
of nuclear equation of state (values of incompressibility and density dependent behavior of the 
nuclear symmetry energy ) and the lepton fraction. Nonetheless, the main feature is unaltered, 
especially up to low values of n [27]. Moreover, the contribution on the lepton energy, as is 
presented in Fig. 8, originates mainly from electrons while neutrino contribution is smaller (but 
not negligible). 

The present results for the electron and neutrino energy per baryon, can also be accurately 
reproduced, at least for low values of T, by applying the analytical formula used by Onsi et al. 
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[50] where the energy density e; and energy per baryon Ei of the leptons are given by 



In Eqs. f H3|) and fH4|) . C and C are constants. 

The contributions of baryon and leptons on the total pressure are presented in Fig. 9. In con- 
trast to the situation of the internal energy, the nuclear part contribution plays a more important 
role compared with the lepton part. The lepton pressure Pi is comparable to baryon pressure Ph 
up to 72 ~ 0.2 fm~3, but for higher values of n it is significantly small. What is more, the main 
part of Pi originates from electrons compared to neutrinos as presented in Fig. 10. 

As pointed out by Bethe et al. [9], the crucial feature in determining the evaluation of a 
collapsing pre-supernova core is that the entropy per particle is very low, of the order of unity (in 
units of the Boltzmann constant /c^), and nearly constant during all the stages of the collapse up 
to the shock wave formation. Therefore, the collapse is an adiabatic process of a highly ordered 
system. So, since the supernova matter is characterized by a constant entropy and constant lepton 
fraction, we shall also discuss the properties under this condition. This can be done by converting 
the results for isothermal case (T=const) into those for adiabatic case (S'=const) in terms of the 
T — n relation constrained by a constant entropy. 

The T = T(n) relation is constructed by {T,n} values to satisfy ^(n, T)=const in an 5 — n 
diagram. Fig. 11 shows the results for Yi = 0.3 and Yi = 0.4 for 5* = 1. Temperature is an 
increasing function of n. Furthermore, for the same density, the temperature is higher for lower 
values of Yi. The values of T for various values of n are derived, for the two cases, with the 
least-squares fit method and found to take the general form 

T{n) = an!', 

where a = 35.412, b = 0.70481 for Yi = 0.3 and a = 32.35706, b = 0.67694 for Yi = 0.4. The 
results of this study are in very good agreement with those of Takatsuka et al. [27]. The stars at 
lower density denote the {T, n} values for S = 1 and Yi = 0.4 which are derived from Lattimer et 
al. [in]. It is concluded that the temperature increases considerably when moving from the outer 
part of the star to the center in order to maintain a constant value of the entropy per baryon. 

By applying the relation T — n presented in Fig. 11, the fractions Yi of isothermal case (see 
Fig. 3) is converted into the isentropic one for S = 1. The population of constituents is plotted 
in Fig. 12 as a function of n for = 0.3 and = 0.4. The most striking feature of the results is 
the slight dependence of the fraction Yi from the baryon density n (the same behavior and similar 
results have been found in Ref. [27j). To sum up, during the adiabatic collapse of a supernova, the 
population of the constituents (neutrons, protons, electrons and neutrinos) are almost the same 
independent of the density n. 

The entropy contributions (for 5* = 1) from the constituents are presented in Fig. 13. The 
contributions of the baryon are increasing functions of n, the contribution of electrons is increasing 
function of n, while 5*,^^ is almost independent of density. Roughly, neutrons, protons, electrons 
and neutrinos contribute to S (= 1) by about 50, 30, 17, 3% (for Yi = 0.3) and 45, 30, 20, 5% (for 
Yi = 0.4) respectively. Our results are quite consistent with those of Ref. [27] . 

Fig. 14, displays the internal energies per baryon of respective components Ei versus n for 
Yi = 0.3, Yi = 0.4 and S = 1. The main conclusions of the results are similar to those of 
the isothermal case (see Fig. 7). The lepton energy (which mainly originated from electrons 
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contributions) dominates in the internal energy of the supernova matter even for high values of 
the density n. The nuclear contribution on the internal energy dominates only in high values of n 
(depending on the lepton fraction Yi). 

Finally, in Fig. 15 we compare the EOS's between supernova matter and cold neutron star 
matter. The case for supernova matter corresponds to = 1 and = 0.3. It is thus clear that 
the internal energy Et^t of supernova matter (SM) is remarkably larger than that of neutron star 
matter (NS). As far as the nucleon part is concerned, the in SM is slightly lower than 
that in NS due to the large energy gain in symmetry energy (see also |22j)- However, the lepton 
contribution on the internal energy Ei is remarkably larger in SN matter compared to NS matter 
due to the effect of a large lepton fraction, that is, a large kinetic energy of abundant leptons. High 
temperature also contributes to the stiffening, but it is less effective than the high lepton fractions 
(see also Fig. 7). The present results also correspond well with those presented by Takatsuka et 
al. [271 EH] a few years ago 

4 Summary 

The evaluation of the equation of state of hot nuclear matter is a major challenge for nuclear physics 
and astrophysics. EOS is the basic ingredient necessary for studying the supernova explosion as 
well as for determining the properties of hot neutron stars. The motive for the present work has 
been to apply a momentum-dependent interaction model for the study of the hot nuclear matter 
EOS under /^-equilibrium. We have calculated the lepton fractions by applying the constraints for 
chemical equilibrium and charge neutrality. The internal energy and also the pressure have been 
calculated as functions of baryon density and for various values of temperature. Special attention 
has been dedicated to the study of the contribution of the components of /3-stable nuclear matter 
on the entropy per particle, a quantity of great interest in the study of structure and collapse 
of supernova. We have presented and analyzed the contribution of each component. Finally, we 
have presented the EOS of /5-stable hot nuclear matter, by taking into account and analyzing the 
contributions to the total pressure of each component. The above EOS can be applied to the 
evaluation of the gross properties of hot neutron stars i.e. mass and radius. 
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Figure 1: The difference S{n,T,I) — S{n,T,I = 0) as a function of P at temperature T = 10 
MeV and T = 30 MeV, for two baryon densities. 
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Figure 2: a) The contribution of the various terms V^, V^, and the total potential energy 
density Vint as a function of the baryon density b) The momentum dependent term V'~^ as a 
function of the baryon density at temperature T = 0, 10, 30 MeV. 




Figure 3: Electron Ye and neutrino Y^^ fraction versus baryon density for various values of T for 
total lepton fraction a) Yi = 0.3 and b) Yj = 0.4. 
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Figure 4: a) The nuclear symmetry energy for three different parametrization (equation fllip ) of 
the interaction part with the results of reference [H] (see text for more details) and b) the proton 
fraction Yp versus baryon density for cold neutron star matter (down curves) and supernova matter 
(up curves) for the three different parametrization of the nuclear symmetry energy. 




Figure 5: Electron /ie and neutrino /i^^ chemical potentials versus baryon density for various values 
of T for total lepton fraction a) = 0.3 and b) Yi = 0.4. 
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Figure 6: Contribution to the total entropy per particle of a) baryons (Sb), b) leptons (Si) and c) 
the total entropy (Stot) versus the baryon density for various values of T for total lepton fraction 
Yi = 0.3 and Yi = 0.4. 



18 




Figure 7: Contribution to the total energy per particle of baryons (Eh), leptons {Ei) and the total 
energy (Etot) versus the baryon density for various values of T for total lepton fraction a) Yj = 0.3 
and b) y; = 0.4. 
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Figure 8: Energies per particle of respective components Ei versus the baryon density n for the 
specific case with T = 30 MeV and Yi = 0.4. 
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Figure 9: The pressure of baryons (Pb), leptons (P;) and the total pressure (Ptot) versus the baryon 
density for various values of T for the cases a) Yi = 0.3 and b) Yi = 0.4. 




Figure 10: Pressures of respective components Pi versus the baryon density n for the specific case 
with T = 30 MeV and Yi = 0.4. 
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Figure 11: Temperature T-density n relation with Yi = 0.3 (solid line) and Yi = 0.4 (dashed line) 
for S = 1. Stars denote the results for the case with S = 1 and Yi = 0.4, extracted from the 
results by Lattimer et al. 
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Figure 12: Fractions of respective components Yi as functions of the density n with 5* = 1, for 
Yi = 0.3 (solid lines) and Yi = 0.4 (dashed lines). Stars denote the results of Yp for the case with 
5* = 1 and Yi = 0.4, extracted from the results by Lattimer et al. [19]. 
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Figure 13: Entropies per baryon of respective components Si contributing to the total entropy 
5" = 1 of supernova matter with total lepton fraction &) Yi — 0.3 and h) Yi — 0.4, as functions of 
the density n. 




Figure 14: Internal energies per baryon of respective components versus n with S — 1 for the 
cases &)Yi = 0.3 and b) Yi = 0.4. 



22 



250 I — I — I — I — I — I — I — I — \ — I — I — I — r 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 

n(fm"^) 



Figure 15: Internal energy per baryon versus density n for dense supernova matter (SM) in 
comparison with that of cold neutron star matter (NS) by applying the same model. The case of 
supernova matter corresponds to 5" = 1 and Yi — 0.3. The contribution of each species is plotted 
separately. 
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